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TECHNICAL NOTE 3157 


METHOD FOR CALCULATION OF COMPRESSIBLE LAMINAR 
BOUNDARY LAYER WITH AXIAL PRESSURE 
GRADIENT AND HEAT TRANSFER 


By Paul A. Libby and Morris Morduchow 
SUMMARY 


A rapid and sufficiently accurate method, for most practical purposes, 
of determining laminar -boundary-layer characteristics in flow with a given 
free-stream Mach number and given velocity distribution at the edge of the 
boundary layer is presented. The method can be easily applied to flow with 
zero pressure gradient for any (constant) Prandtl number of the order of 
unity and any given temperature distribution along the wall. Numerical 
examples are given to illustrate the method and the satisfactory accuracy 
obtained. For flow in an axial pressure gradient, the method can be 
applied for a Prandtl number of unity and any given uniform wall temper- 
ature. The methods developed here are based on an application of the 

, integral method to both the momentum and energy equations, in 
conjunction with a sixth-degree velocity profile and a seventh-degree 
stagnation-enthalpy profile. A single boundary-layer thickness and one 
of the coefficients in the thermal profile are the parameters in this 
two-parameter method. 


INTRODUCTION 


The aim of this report is to present a relatively simple method, 
sufficiently accurate for most practical purposes, of calculating the 
laminar-boundary-layer characteristics in the compressible flow over а 
given object with heat transfer at the wall. 


The method is based on the extension of the Kírmén-Pohlhausen method 
to sixth-degree velocity profiles and seventh-degree stagnation-enthalpy 
profiles. The use of sixth-degree velocity profiles is in accordance with 
the conclusions of reference 1, wherein it was found that such profiles 
can usually be expected to lead to results of adequate accuracy without 
much increase in computational work. Such profiles have been applied 
with satisfactory results for compressible flow over a flat plate with 
heat transfer (ref. 2) and for compressible flow in an axial pressure 
gradient without heat transfer (refs. 5 and kh). 
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The present analysis consists of two main parts. In the first a 
simple approximate general solution is developed for flow without axial 
pressure gradient, with any given Mach number (neglecting hypersonic y 
effects) and given distribution of temperature along the wall. The 
Prandtl number is also left arbitrary, although it is required to be of 
the order of magnitude of unity. The coefficients of specific heat, as 
well as the Prandtl number, are assumed as constant, while the viscosity 
coefficient is assumed to vary linearly with the temperature, the propor- 
tionality factor being chosen in accordance with the Sutherland relation, 
ав in references 1, 2, 5, 5, and 6. Іп the second part of this analysis, 
а simple approximate general solution is developed for flow in a pressure 
gradient with any given constant wall temperature, free-stream Mach num- 
ber, and velocity distribution outside the boundary layer. The Prandtl 
number is, however, assumed to be unity. 


A brief indication of the pertinent literature may be worth while. 
The flow without axial pressure gradient (e.g., flow over a flat plate) 
has, of course, already been treated to a considerable extent. Kármán 
and Tsien (ref. 7) derived solutions for a Prandtl number of 1, constant 
wall temperature, and viscosity coefficient proportional to a power of 
the temperature. Chapman and Rubesin (ref. 5) have obtained exact (numer- 
ical) solutions for variable wall temperature and a linear viscosity- - 
temperature relation. Van Driest (ref. 8) has accurately calculated a | 
series of curves showing skin friction, heat-transfer coefficients, and 
velocity and temperature profiles for a constant wall temperature, using v 
the Sutherland viscosity-temperature relation. Klunker and McLean 
(ref. 9) have calculated the boundary-layer characteristics using the 
actual variation of Prandtl number, specific heat, and viscosity coeffi- 
cient with temperature, while Young and Janssen (ref. 10) and Moore 
(ref. 11) have done likewise with the aid of a differential analyzer, 
Moore taking dissociation into account. 


In each of the above references, the calculations, though exact, are 
fairly tedious, while, except for reference 5, the wall temperature is 
assumed as constant. The advantage of the method presented in the present 
report is that the calculations can be performed with relatively little 
difficulty for a prescribed wall-temperature distribution, Mach number, 
and Prandtl number. Although the results thus obtained will not be exact, 
they will usually be sufficiently accurate in practice. 


In contrast with the literature on flow over a flat plate, the liter- 
ature on compressible boundary-layer flow in a pressure gradient with heat 
transfer still appears to be limited. А good survey of work done on 
incompressible flow can be found in reference 12. Kalikhman (ref. 15) w 
has treated the compressible case by applying the Kármán-Pohlhausen method M 
to both the momentum and energy equations. Ferrari (ref. 14) has recently | 
presented an accurate, but elaborate, method of calculating the boundary- v 
layer characteristics for arbitrary Prandtl number and uniform wall 
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temperature. The present analysis has, once again, the advantage of 
presenting а general approximate solution which involves relatively little 
computational difficulty. However, the Prandtl number is here restricted 
to unity. 


The method of analysis used in this investigation merits some comment. 
As in references 15 and 15, stagnation-enthalpy, rather than temperature, 
profiles are used. The use of stagnation-enthalpy profiles leads to 
significant mathematical simplifications, while it can be shown, more- 
over, that for a Prandtl number of 1 the equations used here automatically 
lead to what are known to be exact integrals of the energy partial- 
differential equation. As has already been stated, the velocity profiles 
and stagnation-enthalpy profiles are respectively of sixth and seventh 
degree here. Іп reference 15, both profiles are of the more customary, 
but usually less accurate, fourth-degree type. Finally, it is of interest 
to observe that in the present analysis only a single boundary-layer 
thickness is used, instesd of & dynamical and & thermal boundary-layer 
thickness, such as used in references 15 and 15. This, however, does not 
necessarily impose any undue restrictions on the thermal profiles, since 
the latter have here been permitted to contain an additional coefficient, 
not determined in advance by the boundary conditions. This additional 
coefficient, which must be determined by the differential equations, 
replaces the thermal boundary-layer thickness as & parameter. This 
explains why the stagnation-enthalpy profiles are here of one degree 
higher than the velocity profiles. For Prandtl numbers near unity, as 
is the case for air, it has been found (cf., for example, refs. 13, 15, 
and 16) that the ratio of dynamical to thermal boundary-layer thickness 
is close to unity. Consequently, the present analysis should yield 
sufficiently accurate results for such cases. The advantages of using 
& single boundary-layer thickness are that the equations thus obtained 
are somewhat simpler and that 1% appears somewhat easier to solve for 
the additional coefficient than for the second boundary-layer thickness, 
especially for flow without a pressure gradient.l 


This work, carried out at the Polytechnic Institute of Brooklyn 
Aeronautical Laboratories, маз sponsored by and conducted with the 
financial assistance of the National Advisory Committee for Aeronsutics. 
The authors are happy to acknowledge the aid of Mr. Richard Grape in 
carrying out the numerical calculations and in checking the equations. 


This is due primarily to the fact that in the two-thichness method, 
the ratio of the thicknesses appears in the equations in & complicated 
manner. 
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constants in equations (50) and (51) р 
constants defined by equations (63) 

coefficient of т] іп velocity profile (eq. (20)) 

constant "average" value of ас 

coefficient of т] in stagnation-enthalpy profile (eq. (21)) 


constant "average" value of bj 


proportionality factor in temperature-viscosity relation 
(eqs. (9) and (10)) 


constant "average" value of С. 


i ди /д 
Л (udu/dy)}, ах 


average skin-friction coefficient, i 
2 
(1/2) p, д 


local skin-friction coefficient, (20/35) /(рым„2)2 


coefficient of ed in given wall-temperature distribution 
(eq. (A1)) 


specific heats at constant pressure and at constant volume, 
respectively 


coefficient of Ё іп equation (А5) 
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Fy, Fo, Fs integrals appearing in equations (11) and (12) and defined 
thereafter (cf. also eqs. (61)) 

P. E, constant "average" values of F4 and Р. 

ғ, ‘functions of т appearing in equation (А5) and defined there 


б1(&) = Ho/83 (8); for Np. = 1, G = To/Te (cf. also eq. (46)) 


Gos Gz (),G, (8),6, (8),6 (Е) parameters defined in equations (33), (38b), 
(45b), and (58) 


H stagnation enthalpy, (42/2) + cT 
He stagnation enthalpy corresponding to equilibrium well 
temperature 
h local heat-transfer coefficient (eq. (40)) 
k coefficient of heat conductivity 
L characteristic length 
M Mach number 
ш,3,2 constants defined by equations (67a) 
N Nusselt number 
Nu 
Nop Prandtl number, Hep К 
р Static pressure 


heat-transfer rate, -(k oT/dy) o 


R gas constant, Cp 7 су 

Ry, Reynolds number based on 1, рауы, 

к, Reynolds number based on x, ру, 

S Sutherland constant; 216° R for air (cf. eq. (10)) 


T absolute temperature: 


T 


^ = Ry (84/1) ° 


Ц 


у 


Subscripts: 


О 


1 
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equilibrium wall temperature for zero heat transfer 


transformation variable, defined by equation (13) or (14) 


velocity components in x- and y-directions, respectively * 
coordinates par&ilel and normal to surface, respectively 
constant parameters, depending on М. and defined after 

equation (Ak) 
ratio of specific heats, Cy [Cy y = 1.4 for air 
boundary-layer thicknesses іп xy- and xt-planes, respec- 

tively (cf. also eqs. (15) and (16)) 
recovery factor, defined by equation (42) 
coefficient of viscosity 
kinematic viscosity, u/p у 
dimensionless distance along wall, x/L ‚ 


mass density 


dimensionless variable, %/5, 


values at wall; for example, То 
local values at outer edge of boundary layer; for example, 
Ту апа М Е | 


values in undisturbed free stream; for example, М, (Also, 
ef. footnote 2) 
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GENERAL EQUATIONS 
The following equations describe the steady, two-dimensional, laminar- 
boundary-layer flow of a compressible gas along а slightly curved wall: 


Momentum equation in the x-direction: 
рии, + руш, = -Py + (uuy) (1) 
у 


Momentum equation in the y-direction: 


= О 2 
Py (2) 
Continuity equation: 
(ou). + (оу), = 0 (3) 
Equation of state: 
p = pRT (4) 
Energy equation: 
pucp Ty T руср Ту = ир, t (erg) + ищу” (5) 


Equation (2) implies that the static pressure within the boundary layer 
is the same as the static pressure in the potential flow just outside the 
boundary layer. Thus, for slightly curved walls 


-р, = -ар/ах = руч (au, /àx) (6) 


where the subscript 1 denotes local conditions at the outer edge of the 
boundary layer. It may be noted that the assumption that the pressure 
in the boundary layer is equal to the potential flow pressure must be 
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modified at high Mach numbers and very low Reynolds numbers (cf. ref. 17). 
The methods to be used in the analysis of such flows are as yet not clear. 
As is usual in aeronautical problems the coefficients of viscosity and 
heat conductivity are considered to be known functions of the temper- 
ature, while the coefficient of specific heat at constant pressure is 
assumed constant. 


It is convenient in the method of analysis to be presented here to 
rewrite equation (5). If equation (1) is multiplied by и and added to 
equation (5), it is readily possible to put the resulting equation in 
the form 


tent, on) = fafi - Gn) бе т 


У 


provided the Prandtl number is assumed constant. For most gases the 
actual change of Prandtl number with temperature is sufficiently small 

to justify this assumption, especially since the experimental determi- 
nation of the heat conductivity appearing in the Prandtl number is subject 
to considerable error (cf. ref. 18). 


Before discussing a method of treatment of equations (1), (3), (4), 
(6), and (7) the mass density and coefficient of viscosity will be written 
explicitly in terms of the temperature. From equations (2) and (4) one 
obtains the equation 


p/P; = T. m (8) 


For mathematical convenience the viscosity~temperature relation suggested 
in references 5 and 19 is extended here to the case of pressure gradient. 
Thus, within the boundary layer 


ry Bas ‘= с(т/т,) (9) 


\ 


where 


VES (т./т„)/ : (т, + 5) / (v, + S) (10) 


With this choice of C the approximation of equation (9) to the actual 
viscosity-temperature relation is good in the neighborhood of the wall, 


2Ү 
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where the skin friction and heat-transfer properties are determined, and 
is, for most problems, considered satisfactory throughout the boundary 
layer. It may be pointed out that this viscosity-temperature relation 
is equivalent to using u/ui = c (m/m) where C is similarly defined 


in terms of То and T}, if it is assumed that the ratio р/н, is 
calculated by means of the exact Sutherland relation. 


By integrating equations (1) and (7) with respect to +t over the 
boundary-layer thickness t = O to + = 5+ and using the boundary con- 


ditions u-v-0O at t= О, together with smooth transition of the 
velocity and temperature profiles to their local main-stream values, the 
following differential-integral equations are obtained: 


(5, B) + A4ES + Еу (loge 01)' + (log, ш)" ғ, + d " L— ња)», е 
се) Бера) (8/7) ТЕЙ! an 
(2/2) AP EA FI е, бов, 4); + (Log, 1l = 
(fe) (a/e) 6 Gne [tr aa 
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= 5/5, 


The prime denotes differentiation with respect to E, and % is the 
Dorodnitzyn variable defined by the transformation 


J 
t = T. /T 15 
J, (4/9 % as) 
or by the inverse 


y -Г (т/т,) at (14) 


Moreover, 


54 = [А (т, /2) dy (15) 


or 


5 JU (т/һ,) at (16) 


The quantities 0, /P and M, in equations (11) and (12) are 
со 


related to ча» which 15 a function of Ё prescribed by the potential 


flow about the body in question. Thus, in accordance with the usual 
isentropic flow relations, 


?In applying these relations and interpreting infinity conditions, 
the presence, in the supersonic case, of attached or detached bow waves 
should be considered. "Infinity" conditions should then be taken as con- 
ditions just behind such shocks on the streamline which in potential flow 
would form the body under consideration. 


NACA TN 3157 


yu 1/(y-1) 


оу Po = (т, /т) = 41 + (у - 1) (м 2/2) Е z (% Іш) і (17) 


1/2 


му = (ufa) le  -D ЖЕР : md | (28) 


Approximate solutions to equations (11) and (12) can be found by 
assuming u/u, and H/H, as functions of т which satisfy certain 


boundary conditions at т = O. and т = 1 to be discussed in detail 
later. In these profiles one parameter in addition to A is left unde- 
termined after the boundary conditions are satisfied, and these two param- 
eters are considered as the two unknown functions of t€ to be determined 
by equations (11) and (12). 


ln setting up the final differential equations, examination of equa- 
tions (11) and (12) indicates that, in addition to selecting the profiles 
for ола and H/E with their appropriate boundary conditions, one 


must prescribe (21/2. 6), Мо and Npp. Moreover, the values of the 


со? 
fluid properties at infinity must also be specified. From the solutions 
to equations (11) and (12) all quantities in the &t-plane can easily be 
found; the transformation, equation (14), can be applied to obtain the 
corresponding quantities in the xy or physical plane. For this purpose 
equation (14) may be conveniently written as 


y/L = (54 L) E (ҥ/н,) [ + (y - 1) EYE] = (y = 1) (9/2 (а/ш) dt 


(19) 


The boundary conditions that the velocity and stegnation-enthalpy 
profiles are to satisfy will now be discussed. As indicated in the 
"Introduction" a sixth-degree polynomial in т for the velocity and a 
seventh-degree polynomial for the stagnation enthalpy are suggested on 
the basis of the critical study of the use of integral methods in 
compressible-boundary-layer analysis presented in reference 1. Thus, 
let 


12 


All of the а, and 5, 


boundary conditions at 
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6 
= ат 20 
ш/а] 2 j (20) 
T 
д 
H => br (21 
f^. = 24 ^s | 


coefficients except one are determined from 


т = О and т-1. In setting up the differential- 


integral equations (11) and (12), the boundary conditions which must be 
specified to obtain a unique and exact solution to the governing partial- 


differential equations 
At т = 0, 


and аб T = l, 


have been applied. Іп the ёт-ріапе these are: 


ufu = v/a = О 
(22) 
Н/Е = E5 (&) 
ufu = H 1 = 1 
(23) 


The boundary conditions specified here are the usual ones іп heat-transfer 


problems. It is possible to specify instead of он), either ko 


1 о 


or some relation between these two quantities as in the case of radiation 


cooling. 
is left arbitrary and 


In principle there is no loss in generality here, since Ho / Hy 


nothing is specified with respect to 8/83) | : 
T 
о 


4 


4: 
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In addition to requiring that the assumed profiles satisfy boundary 
conditions (22) and (23), it is customary in the integral methods (cf., 
e.g., ref. 1) to select the profiles so that at the wall and at the edge 
of the boundary layer some of their derivatives have the same values as 
those which an exact solution to the original partial-differential equa- 
tions under boundary conditions (22) and (23) would yield. For the pro- 
files of equations (20) and (21) these additional conditions can be 
obtained from equations (1) and (7) and from differentiation of these 
equations with respect to +t. Іп this way one obtains the following 
boundary conditions to supplement those of equations (22) and (23): 

At т = 0, 


обл.) (Уа), = Mer.) sf) + © - al а)" 9 


CO) quA AR (25) 


2 


1 : x, ) (у - i ler), | (26) 


bre nis] өы), 
эу» Pu OP) E)" = “(АА n)... 
М-ы) xs) ef), о en 


At т- 1, 


(71) 7 (Ун) ттт І (У "m. І Uu M a = 
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With equations (22) to (28) fourteen of the fifteen coefficients in 
the velocity and stagnation-enthalpy profiles of equations (20) and (21) x 
can be determined in terms of A, of the fifteenth coefficient (to be 


chosen as bj), of the prescribed functions (u fù.) (£) and (8/81) 6), 
and of the prescribed parameters M, and Np. 


As explained in the "Introduction," it will be observed that only 
one boundary-layer thickness 68, is here assumed and that thus no 


distinction is made between a velocity and a temperature boundary-layer 
thickness. This, however, does not impose an undue restriction on the 
shape of the temperature profiles, since the latter (see eq. (21)) are 
here permitted to retain an additional parameter to be determined, not 
by any boundary conditions, but by the differential equations (viz, 
eqs. (11) and (12)). 


Examination of equations (24) to (28) indicates that considerable 
simplification is achieved in the boundary conditions, and therefore in 
the analysis, under two special cases: (a) The flow without axial pres- 
sure gradient (ux fu, = 1), which also corresponds to the practically 


interesting E of the supersonic flow over a thin wedge, as well as 

to the subsonic (or supersonic) flow over a flat plate and (b) the air- 
foil with pressure gradient, and constant wall temperature, in а fluid 
with Prandtl number equal to unity. Although the assumption of constant 
wall temperature does not appear to be exactly realized on practical 
wings, it appears desirable to investigate the boundary-layer character- Б 
istics under these assumptions since the pressure-gradient, heat-transfer 

case in a compressible flow has not yet been analyzed entirely satisfacto- 

rily from the standpoint of simultaneous accuracy and simplicity. 


FLOW WITHOUT AN AXIAL PRESSURE GRADIENT 


In this section a general solution based on the previously derived 
integral-differential equations will be obtained for the special case of 
flow without axial pressure gradient. The accuracy of this solution is 
then investigated by comparing it with more exact solutions. 


General Case 
In the case of the flow without axial pressure gradient and with : 


general but constant Prandtl number and variable wall temperature the 
integral-differential equations (11) and (12) become: 
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О 


(5/8) + Еул = С (4), (29) 


(F5/2) At + FUN (с/в) (е) | (50) 


Неге M = №, Ту = Ty; and so forth. The subscript о will be used 


to denote the free-stream conditions for this case. 


It is convenient to consider <A and by as the unknowns to be 


determined by equations (29) and (30). ‘Then the velocity profile satis- 
fying boundary conditions (22) to (28) is 


6 


ufu, = 2r - srt + бт? ~ Or (51) 


While the b j coefficients of the stagnation-enthalpy profile become: 


bo = Н/Е, = Gy 

bo = 2G, 

55 = 205 

by = 25(1 - 61) - 2062 - 86, - 20b, (52) 
bs = -8} (1 - Cy) + НОбо + 12@5^ + h5bi 


56 = ТО(1 ~ Gy) - 5062 - 885^ - 5601 


b; = -20(1 - 6) + 862 + 205^ + 10%] 
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where the G's are known functions of € defined by: 


Cy = (E/E) (е) 
бр = (1 - Np.) (7 - an УҺ + (у - 1)(м„/ 2) (55) 
G5 = (1/6)Np. оу С 


From the definition of Ру and from equation (51) it ів found that 


Ру = 985/9,009 (5%) 
Similarly, 


РУ = (51/126) (1 5 вт) - (955/180,180)83^ - 
(302/9,009)Gp ~ (821/12,012)b4 (35) 


With the value of Е] given by equation (34), ACE) can be found 
directly from equation (29). ‘Thus with the requirement that A(O) = O, 


A = hey Ary (36) 


where” 


S 
а-/ С dé 
О 


Equation (30), which gives %4(Е), becomes: 


Ёз! - (у н.) (с/с o3 + (ғ /2) (с/с) = © (37) 


ITt may be noted that in this analysis the temperature-viscosity 
proportionality factor С, unlike in reference 5, can be permitted to 
vary along the plate (1.е., С = C(E)) without any mathematical difficulty. 


ЗҮ 
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^ The solution of equation (37) satisfying the requirement that b be 
finite at E = O is | 


- pë RC 
bj = Cy E c, ~G, (8) at (38a) 


where 


B1 = (1/2) + (985/2,463) (1/%ь,) 


G, (Е) = (12,012/821) (буе 1-6,)(c/c,} - 
Л (1 - 9) (6/8) TT 


э вул) + (955НЬ, /29›550) (б/г " (осу с) _ 


К (955Npr /29,550) (су /с)су" - (151/9,009) (6/63) i 


From equation (38a) with any prescribed temperature distribution 
G (£) at the wall, bi = bi(t) can be readily determined. From equa- 


tions (32) the remaining Ы; coefficients can then be found as functions 
of t and of the prescribed flow parameters, for example, М>, and М. 
With the а; апд b; coefficients and A all determined, the velocity 


and stagnation-enthalpy distributions throughout the €t-plane аге known. 
By application of equation (19), these distributions in the xy or phys- 
ical plane are easily determined. 


The average skin-friction coefficient С. for the length 1 can be 
calculated from equations (31) and (36) ав 


L 
where Cy -/ CGE. 
| О 
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The local heat-transfer coefficient h is usually defined as 


hz -(kmy) | Дз - Те) (40) 


Using solution (36) for A, equation (40) yields: 
h = -0.165 үнү c/o - Gy 


E 1 + (у - 1) ( м, eel + Я - 1)( 718 (41) 
| where 1 is the recovery factor defined by the equation 
T е/ Тү = 1 + (у - 1)( (м, 2) = (8/81) 1 + (у - (м2 В) (42) 


The (nondimensional) Nusselt number Ny, is thus 


0.165 c Af, ) Ву, Eby 
E E + (y -= SOWAN z (y - 42) 


The value of тү .follows by setting by = О and solving for the 
constant wall temperature, as determined by бу. Equations (33), (38), 
and (42) with у = 1.4 are thus found to imply: 


(43) 


Nyy = х/к = - 
G 


n=l- 0.272(1 - Қы.) (44) 


NACA TN 3157 19 


It should be pointed out that equation (44) yields for Np, = 0.72 


а value of т = 0.924, which must be compared with the value of 0.845 
given by the mathematically more exact method of Chapman and Rubesin. 
Although this result must be considered a consequence of the approximate 
method of analysis used here, calculations involving rates of heat trans- 
fer (such es Nusselt number) for given values of the wall temperature 
will not be directly affected by the equilibrium wall temperature and, 
therefore, will not be appreciably affected by any inaccuracy in 1. 

This is borne out by the numerical example and the comparison with an 
exact solution, given in a subsequent section here, with NS. = 0.72. 


It may nevertheless be convenient for certain applications to express 
the wall temperature in terms of the predicted equilibrium temperature. 
With the additional simplification of replacing variable C by an average 
value C equation (43) becomes: 


mJ МА, = 0.297 VE h - б/т.) кш J i G5 dé (45a) 


where 


В1-1 


G5 = 5 Е - (То те) - а + (бо, 055/152, 615)» (To/TeJ § T 
(10,026/152,615) М>. (Tof Te) nee (455) 


In deriving equations (45a) and (45b), the following relation between бу 
and тт (cf. eq. (h2)) was used: 


ву = (To/Te) R + (y - 1) м,2 f^ : + (y - 1) (м2 2) (46) 


"Because of the inaccuracy in т, it will be found, of course, that 
when the wall temperature has the same value as the true equilibrium 
temperature, the present equations will not lead exactly to a zero heat 
transfer at the wall. However, in such a case, it will be found tat the 
calculated rate of heat transfer will be relatively smail. 
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It should be observed bag the recovery factor т does not appear in 
solution (45a), although а. is prescribed. 


If it is desired to calculate the actual rate of heat transfer (4 
(in British thermal units per second per square foot) for a given 
ratio T [ты from equation (45a), then the Nusselt number as determined 
by equation (45a) should be multiplied by (k,/X T, - Tej» where T, ің 


the actual equilibrium wall temperature (as given, e.g., by exact theory 
or by experiment). Іп practice, however, the wall temperature To (or 
the ratio Т./Т,)» rather than the ratio of wall temperature to equilib- 
rium temperature, will usually be prescribed. In that case, the rate of 
heat transfer can be directly calculated (without any intermediate use 
at all of the recovery factor т) from the expression: 


а = -(k òy), = -C A (1 т 2-2 m’) T ii т) 


The quantity bq can be calculated from equations (58а) апа (580) for 
flow over a flat plate, with the observation that 


a= (ыл) неу ee) 


Equation (47) has been derived with the use of the following relation 
between temperature and stagnation-enthalpy profiles, resulting from the 
definition of the stagnation enthalpy Н: 


iuum) teg) бә 


In order to determine the velocity and temperature profiles in the 
physical (xy) plane, the inverse transformation given by equation (19) 
must be used to find the coordinate у associated with the variables Ё 
and T, while the x-coordinate may be found directly from the definition 


of E. 
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If the wall-temperature distribution is prescribed as a polynomial 
in the distance € along the wall, then the Nusselt number and temper- 
ature profiles, following from equations (45a), (45b), and (49), can be 
directly calculated, without any quadratures, from the explicit expres- 
sions given in the appendix. 


It is of interest to observe? the wall-temperature distribution 
required for zero heat transfer according to equations (45a) and (№5Ъ). 
Thus, by setting G5 = 0, а second-order linear differential equation 


in То/Те is obtained, whose general solution is: 


-1/2 -T.65/N 
2 РТ + АЕ + ВЕ ру (50) 
е 


H 


where A and B are arbitrary constants. Equation (50) implies that 
zero heat-transfer conditions will be satisfied not only by a uniform 
wall temperature Te but also by а wall temperature varying inversely 


as а certain power of the distance Е along the wee In particular, 
-1/2 


the heat-transfer rate vanishes for 1 - 7) с Ё (regardless of 


the Prandtl number). This result was also obtained by Levy (ref. 20) for 
incompressible flow, using an analysis restricted to power-law wall- 
temperature distributions. The wall temperature according to this distri- 
bution would be infinite at Ё = О, but Eckert (ref. 21) has reasoned 
that a heat-line source or sink placed at the region where Ё = O might 
result in ә temperature distribution of this type. ‘The case of upstream 
cooling such as analyzed in reference 22 is the physical counterpart of 
the abstraction of this line source or sink. Equation (50) indicates 
that there is another type of wall-temperature distribution for zero heat 
-7 -65/Npr 
transfer, namely, 1 - (T,/T,) = & . This is not included in 
reference 20, where negative powers of & for T (£) lower in absolute 
value than 7-63/N,.. (approximately 10 for air) were assumed. It should 
be noted that the results obtained in the present analysis are valid for 
compressible flow (i.e., for any Mach number). Finally, it may be noted 
that temperature distributions of the form of equation (50) are not 
included in polynomial distributions, such as those to which reference 5 
is restricted. 


ӘТһе observations in this paragraph were pointed out to the authors 
by Dr. M. Bloom. 
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Constant Wall Temperature and Np, = 1 


A special simple case of the above solution for flow over a flat 
plate is of interest here. This is the case of constant wall temperature 
and Prandtl number equal to unity. This case will now be shown to give 
а particularly simple check on the accuracy of the method presented here. 


It is well known that if the Prandtl number is assumed equal to unity 
then equations (1) and (T) imply for the flat-plate case that 


H= А + № (51) 


where A and B are arbitrary constants. Equation (51) is an exact 
stegnation-enthalpy relation which satisfies the boundary conditions of 
constant wall temperature both with and without heat transfer. Therefore, 
one test of the approximate method presented above is its automatic 
reduction to this exact relation, if the conditions under which this 
relation is known to be exact are imposed. This test will now be applied. 


For the case of constant wall temperature and & Prandtl number of 
unity, equations (58Ъ) and (55) imply: 
Gy = (4,433/2,463) (1. - G,) (1/8) 
Hence, equation (38a) yields: 


b, = 2(1 E G4) (52) 


Furthermore, from equations (32), 


bo = bz = О 
b, = -5(1 Е с.) 
55 = 6(1 - Gy) (53) 


bg = -2(1 x с 


be = О 
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When these bj coefficients are substituted into the stagnation-enthalpy 
profile, one obtains 


H/I} = бу + (1 = ст) (2r - Sr + 67? - 2т6 ) (54) 


end thus, by comparison with equation (31), 


H/E; = @у + (1 - &) (ufu) (55) 


Equation (55) is exactly equation (51) in nondimensional form, with 


the boundary conditions that H = Но when u- 0 and Н = Ну when 


u=u,- It is thus seen that the reduction to the exact energy equation 


for this special case is automatic and rigorous. It has already been 
shown in reference 1 that the agreement between the results of the inte- 
gral method using а sixth-degree velocity profile in conjunction with 
equation (55) (i.e., for constant wall temperature and Prandtl number 
equal to unity) and the results of the exact solution of reference 5 is 
good. The skin friction and heat transfer are for all practical purposes 
exactly predicted and the laminar-boundary-layer stability limits are in 
good agreement. It can therefore be expected that the more general inte- 
gral method presented here, valid for cases of heat transfer іп а ргев- 
sure gradient, will be accurate in more general cases, since the two 
integral methods are in complete agreement for the special case just 
considered. 


Comparison With the More Exact Method of Chapman and Rubesin 
(Variable Wall Temperature and Np, # 1) 


ln reference 5, Chapman and Rubesin have applied their analysis to 
the case of a variable, prescribed wall temperature. To simplify the 
numerical analysis the proportionality constant C was assumed by them 
to be equal to unity. The Prandtl number was assumed to be 0.72. In 
the notation used here the wall temperature was expressed there as 


т./Т. = 1.25 - 0.858 + 0.5552 (56) 


Thus for equation (Al) (see appendix) со = 1.25, су = -0.85, and 
Со = 0.55. Тһе Nusselt number variation with Ё given by equation (A2) 


is shown in figure 1, along with the results of the more exact Chapman- 
Rubesin calculation. It can be seen that excellent agreement is realized. 
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Tt will be noted from equation (59) with с, -1 that the average skin 
friction obtained by the present equations is in excellent agreement 


with the exact value Ce, / Rr, = 1.328, given in reference 5. 


In сойрагіпе the velocity and temperature profiles predicted by the 
two methods it must be kept in mind that for strict comparison the wall- 
temperature distribution should be precisely the same in the two methods. 
However, since Chapman and Rubesin assumed only the ratio of wall temper- 
&ture to equilibrium temperature and since the two methods do not give 
precisely the same recovery factor, and hence not the same equilibrium 
temperature, it is necessary for this comparison to express the wall 
temperature assumed by Chapman and Rubesin, as well as that used in the 
present report, in terms of the free-stream temperature. 


In figures 2 and 3 the comparisons of the temperature and velocity 
profiles predicted by the two methods are shown for three values of 6 
and for two Mach numbers. It can be seen that excellent agreement is 
realized. 


It may therefore be concluded from these numerical results that the 
momentum-and-energy integral method given here will yield good results 
with little computational difficulty when applied to the case of zero 
axial pressure gradient, variable wall temperature, and general but 
constant Prandtl number of the order of magnitude of unity. 


From а practical point of view, it should be kept in mind that the 
solutions developed here are based on the viscosity-temperature rela- 
tions (9) and (10), which are an approximation to the actual relation 
for air. Because of relations (9) and (10), the results obtained here, 
namely equations (39) and (45a), indicate that, for a fixed wall temper- 
ature, the skin-friction coefficient and the Nusselt number will be 
independent of Mach number. For the Sutherland viscosity-temperature 
relation, however, this will not be quite valid (cf. ref. 8). 


FLOW WITH PRESSURE GRADIENT 


For flow over & curved surface, such as an airfoil, where the local 
velocity distribution  u/u, outside the boundary layer is not constant, 


but may be considered as a given function of the distance x along the 
wall, the ordinary differential equations (11) and (12), unlike the case 
of zero pressure gradient, can no longer be easily solved for general 
Prandtl number and arbitrary distribution of temperature along the wall. 
For such flows, however, considerable mathematical simplifications occur 
when the Prandtl number is unity and the wall temperature is uniform 
(ef. especially eqs. (26) and (27)). 
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The assumption Np, = 1 has the additional theoretical advantage 
that, according Чо equation (44), it automatically leads to the known 
exact value of the recovery factor, namely, т = 1, for this velue of 
the Prandtl number. In fact, the parameter Gy now assumes particular 


physical significance, since for Np, = 1 equation (46) implies: 


бу = Т./т, (57) 


Consequently, the case of Np, = 1 and uniform wall temperature will be 
treated now in detail. It will be shown that for such а case a relatively 


simple approximate solution of the equations can be obtained for any 
given velocity distribution outside the boundary layer. 
General Approximate Solution for 
Npr = i and То = Constant 


For this case, with (Е) arbitrary, the coefficients a 
; ч J 


and b, in equations (20) and (21) for the velocity and stagnation- 
enthalpy profiles can, by virtue of boundary conditions (22) to (28), 
all be expressed in terms of ap and bj, where by remains arbitrary, 


while ap is given by 


2-х 
во = -(1/20) (mm) (ааа + Ев) = e» 58) 


where Ge (E) is a given function. The profiles can then be given in 
terms of ас and b, as follows: 
ufu = (er - 57. + 6т2 - 276) + (ғг/5) Е + 5т2 - 107 + 10r - 376 + 


(6; 6e; Jt + 1072 - рот! + 1577 - 6) (59) 
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H/H = G- + (1 - с.) (9% - ёт? + тоте Е 2017) + 


b с - 20r" + 15т2 - 2618 + 107) (60) 


1 


With expressions (59) and (60), the following explicit expressions 


for Fy; F5; and Fz are obtained: 


кг 
li 


0.1095 + 0.002118, - 0.000622. + 0.000422 (уар fa) а 


2 


0.0000095 (5% fy)? - 0.000155 (зе. fea) 


2 
= 0.595 - 0.500(1 : c4 ) + 0.1075, + 0.02128, - 0.000628," + 


Hj 
№ 
1 


, „ > (6) 
0.0028 (5%; fy) : 0,00015(2 8, fea) Е 0.0000095( уг fa) 


„= (2 - 5.) 0-246 - 0.0158, - 0.00181 (ba; ік) : 


»;[o-o665 = 0.003244, à 0.0008 (Ве ік) 


*xj 
| 


With expressions (59), (60), and (61) inserted in equations (12) 
and (12), two ordinary differential equations for М(&) and b,(&) are 


obtained. Although these can be solved numerically for a given ie i DE 


the process may be tedious. A relatively simple general approximate 
solution of these equations will, therefore, be derived. 


Equation (11) can be solved approximately for A by assuming that 
ку and F, can be replaced there by constant "average" values F} 


and Ба over the distance Ё. This is justified by the fact that the 
variable terms there, which are proportional to 2. апа bi; are rela- 
tively small (cf. eqs. (61)). This is equivalent to replacing 82 


and bq by constant "average" values 82 and by for this purpose. 


Al 
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With equation (59) for the velocity profile and equation (58) for а, 


equation (11) then becomes the following linear ordinary differential 
equation in А: 


(FB) + л Ai: fe) + (а /ш) E E i ^ Qe - E) | 
26 (o. р.) (m /т„) бы ш) (62) 


where 94 is & constant, given by: 


Ф =F ВЕ. - 201 - (5. Бо) 


JE 
= А + Db 
(65) 
А = 0.50) + 0.00438 + 0.02528, - 0.001242f 
D = 


0.0905 4 (0.0858 Ж 0.09568) (E, (зов, ) | 


With relations (17) and (18), the solution of equation (62) satisfying 
the condition A = О or finite (if цу = 0 at ё = о) at the leading 


edge (Е = 0) is found to be: 


f зе) E 

Р Cafe) (47%) " 

n (Уғ,)с қ 25 ут. Pi CSE) 
F 7- Fy 
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Equation (64) is similar in form to equations obtained for zero heat 
transfer in references 3 and 23 and for heat transfer, but with fourth- 
degree profiles and two boundary-layer thicknesses, in references 13 
and 15. 


To obtain a general approximate solution for b4 (Е), ^! can first 
be eliminated from equations (11) and (12). Assuming Fy,‘ = ұз = 0 
for this purpose and using equation (58), the following equation is thus 


obtained: 
Е) an ас 
E + fs + (ға (1/5) le (65) 


Noting, in advance, that b, will usually heave a value approximately 
equal to 2(1 - 61), 75,4 can be replaced by 2(1 Š 61) by in the 
relatively small term 0.0000095 (ъа, |б.) appearing in Ру and Fp 

and in the small term 0.0004 1b) “az /бу appearing in Ез. Expressions (61) 
for Fi» Fos 
with coefficients as functions of бу and ал. By substituting these 


and Р. can then be written as linear functions. of by; 


expressions into equation (65), the following quadratic equation in b] 
is obtained: | 


mb? + jb, -1=0 (66) 


where 


E 7 | B. xs 2 5 
m= 222 x 10 les. - h.h96a, + 2.705 28 - 0.2096 22. + 0,01452(=2) + 0.0056 E 
Gy | сү Gy я G1 G1 


0.00052 ag? eoe] 
G2? Gy 


A | | 
j = 0.24602 + 0.0791Ta, - 0.005868а.>2 + = x 10 Еж + 42.260, - 2.06785" + (67a) 


2 | 
0.0465а>? - 7.906 £t 0.9086 Е + 0.067 2) - 0.0252 uc - 0.00284 xd 
| Gy 


= Е е 22 2 
2 eù б) (0.24602 ая). + ol (0.20495 + 0.3G, + 0.021168, - 9) 
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If terms which will ordinarily be negligible are rejected, equations (67a) 
сап be slightly reduced, thus: 


ор, 2 
m = x 10 le.s - В.19бар + = (2.705 - оа) 
Gy G1 
2 = 
j = 0.21602 + ас (0.07917 - 0.005868а..) + = x 10 ы x 


En к (42.26 - 2.06Ta5)as - (T-906 - 0.90868.) 21 (670) 
1 


а. 
Е - с. \ 0.24602 - \ = | - 
1 = 2(1 03) (0.2002 0.01496a,,}< 1 + те 0.10495 + 0.56) + 


а_ (0.02130 - 9.0006216a ) 
2 2 


The solution or equation (66) is: 


: 1 
нея BRE © (68) 


(The physically appropriate root will in general be that which is closer 


to the value 2(1 - Q).) | 


Equations (бі) and (58) represent a simple approximate solution of 
equations (11) and (12) for Np, - 1 and uniform wall temperature 


(& - Constant | . Their application will involve at most numerical 
integration. 


In applying equation (бї), a reasonable "average" constant value 
for аг, for any given (а/%) (Е), С), and Мо, can usually be obtained 


by considering equation (58) for aC /Gi^ and equation (64) for A/C. 
An "average" value by for bj, to be used in evaluating 9, (eq. (63)), 
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can be obtained by considering equation (68) for 9). For objects with 


sharp leading edges, for which A = O at Ё = 0, it will ordinarily be 
found that Ъ- = 2(1~-G,) very roughly. After A(E) has been found 


from equation (64), b,(&) can be easily determined from equation (68). 


With A(t) and Ъ-(Е) determined, the boundary-layer characteristics 
can all be straightforwardly calculated. The local skin friction will be: 


ane e afa jp - bata] (AVR) p tan 


(69) 


The Nusselt number, giving heat-transfer properties at the wall, will be 


Ny, = wee L - (с/ үк) Е e - a) ( (T, /T., ) үш, (70) 


k (то - 


со 


The velocity and temperature profiles follow from equations (59) and (60), 
in conjunction with equation (49) for Т f; and equation (19) for trans- 


forming to the physical plane. -- 


In accordance with conclusions reached in references 1, 2, and 3, 
where the methods applied were essentially the same as in the present 
investigation (viz, sixth-degree profiles, together with а solution of the 
form of eq. (64)), it may be expected that the results obtained by using 
the equations developed in this section will be sufficiently accurate for 
most practical purposes, including stability calculations. However, for 
flow over a blunt-nosed object (i.e., "stagnation flow") or for determi- 
nation of the separation point in an "adverse pressure gradient, these equa- 
tions may have to be modified to yield still greater accuracy (cf. ref. 3, 
where such modifications are shown for zero heat transfer). For stagna- 
tion flow (us fs - KE) in fact, fourth-degree profiles give satisfactory 


б equations (67a) and (67b), the exact expression (58) for a,(€) 
should be used and not any "average" value. 
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accuracy (ref. 5). Consequently, st&gnation flows with heat transfer 
could be calculated by the method of reference 15, where such profiles 
are used for this case. 


Tt should, finally, be observed that equations (64) and (68) have 
been based on the assumption that Fl, F5, and Fz are approximately 


constant, or vary slowly, along the flow. Regarding F4 and Ех, this 


can be interpreted physically as an assumption that the ratio of the 
momentum thickness to the boundary-layer thickness, and the ratio of the 
thermal momentum thickness to the boundary-layer thickness, will be 
approximately constant along the flow in the xt-plane. This approximation 
has been justified by the relatively small variable ал terms (replaced 


here by terms in =>} in expressions (61). In ordinary cases, the &o 


terms will actually be found to be sufficiently small, either individually 
or collectively, so that the approximating assumptions made here will be 
satisfied. In any doubtful case, however, one can in general easily 
check, a posteriori, the validity of the approximate solutions developed 
here by computing the variable ас terms with the approximate solution 


obtained and ascertaining whether these terms have a sufficiently small 
net effect, relative to the other terms, to be replaced by constant 

"average" values. Assumptions corresponding essentially to the use of 
constant "average" values of F4 and Fp have already been applied in 


previous analyses (refs. 3 and 15). 
CONCLUSIONS 


From the analysis presented herein for the laminar compressible 
boundary layer with heat transfer in flow with and without an axial 
pressure gradient, the following conclusions can be drawn: 


1. The boundary-layer characteristics for the flow over a flat plate 
can be easily determined with sufficient accuracy from the equations 
developed here, for a given constant Prandtl number (of the order ОҒ meg- 
nitude of unity) and a given wall-temperature distribution. 


2. For flow with a pressure gradient, the boundary-layer character- 
istics can also be easily determined with sufficient accuracy by the 
equations developed here, provided the Prandtl number is unity and the 
wall temperature is uniform. Here, the velocity distribution outside the 
boundary layer and the free-stream Mach number &re' considered as prescribed. 
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Some modifications of the equations may, for greater accuracy, be neces- 
sary for flow over a blunt-nosed object and for determination of any 


possible separation point. 


Polytechnic Institute of Brooklyn, 
Brooklyn, М. Ү., December 11, 1952. 
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APPENDIX 
POLYNOMIAL WALL-'TEMPERATURE DISTRIBUTION 


If the ratio UN Lu is represented as a polynomial in ё, that 
is if 


d 
Do. e 2 Gui (A1) 
3-5 
then equation (45a) yields: 
M 3 -1 
Хү, М; = em - 2 E (а) (1. = со) -- 


389 (3 + a) d + 25+ 
(60, 059/152, 675) JM. [1 + 15,5%2(4 - У (42) 


With the wall temperature specified as in equation (А1), equation (19) 
together with the use of equations (31), (52), (36), and (49) leads to 
the following explicit expression for у(Е,т): 


(y/2x) В; = V 9, 009€/985 [ -f, + (2/61) + 


(y - 1) м2 Р) (Poe - fz + Baty) + 
И + (у - 1) (м | 2 > e E? fi + В, £54 = 
[ea /(3 + ва E + Bs J T вольт] (A3) 
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where 
f- = т - Tot + 18-2 - тоте T (5/2) 
fo = (4,433/9,852)72|2 - 161? + Зот» - (Ah/T)r2 + =б 
із = - (13/18) ES + 12672 - 19212 + (828/7) 7! Б отт? БА 5018 + 108-17 " 
(1,008/11):8 + 5679 - (72/13)729| 
f, = (22/5) Е - 127° + 20т2 - (90/1) 4- x^ 
f5 = (т* о) о - 32r + hor? - (160/7)т2 + or 


For convenience, the Ғ functions and their first derivatives with respect 
to time, denoted by f,', have been calculated and are plotted in 


figures 4 and 5. 


The temperature profiles in this case can be determined from the 
explicit -expression: 


T/T) = P Es + ға fe) + (у - 1) (м2 Р) (ва! - 25 + 4,” + 
Е + (у - 1)(%2/2)a| > с ја + Byfs'd - 


Es / (3 + ы) |: + Bog + (006/52, 675).°| t (д) 


ay 
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where 


B, = (1/2) + (985/2,463) (1/Npr) 

Bo = E = (1,208/h ‚ь55) (1 = хь.) (2 Pa) 
рур - Npr) 

В, = 28ь,/9Ғі = (6,006/985)8Ь. 

85 = 2 + (20,015/152,675)%ы, 


The velocity profiles пау be found from equations (31) and (А5). 


It may be noted that, if desired, the ratio To f may be preseribed 
аз а polynomial in Ё, thus: 


M 
T/T = 2; et (А5) 


Ав can be seen by comparison of equations (А1) and (А5) in connection 
with equations (46) and (48), the results given by equations (А5) and (AL) 
can still be used, provided that one substitutes 


с, = «s + (у - 1) CVOR (А6) 


there. 
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Figure 1.- Nusselt number distribution Nyu / Rx versus Ё (eq. (42)). 
Flat plate; variable wall temperature; Np, = 0.72. 
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Figure 2.- Comparison of temperature profiles in boundary leyer. 
Variable wall temperature; Np, = 0.72. 


6 


late’ 


(REFERENCE 5) 
velocity profiles in boundary layer. 
; = 0.72. 


© D OR © CHAPMAN & RUBESIN 


NACA TN 5157 


Variable 


ature; Npr 


f 
wall temper 


Figure 5.- Comparison o 


NACA ТМ 5157 


f (т) 


T 


(a) г (т), ғ2(т), апа ғ, (т). 


Figure 4.- Variation of f functions with т. 
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Figure h.- Concluded. 
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Figure 5.- Variation of first derivatives of f functions with т. 
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(b) Ру! (т) and ға (т). 


Figure 5.- Concluded. 


